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ABSTRACT 

We report the detection of the (J,K) = (1,1), (2,2), and (3,3) inversion lines of ammonia (NII3) towards 
the south-western molecular lobe in M82. The relative intensities of the ammonia lines are character- 
ized by a rotational temperature of Trot = 29 ± 5 K which implies an average kinetic temperature of 
Tkin ~ 60 K. A Gaussian decomposition of the observed spectra indicates increasing kinetic temperatures 
towards the nucleus of M82, consistent with recent findings based on CO observations. The observa- 
tions imply a very low NH3 abundance relative to H2, ^(Nlla) « 5 x 10^^". We present evidence for 
a decreasing NH3 abundance towards the central active regions in M82 and interpret this abundance 
gradient in terms of photodissociation of NH3 in PDRs. The low temperature derived here from NII3 
also explains the apparent underabundance of complex molecules like CH3OII and HNCO, which has 
previously been reported. 

Subject headings: galaxies individual (M 82)-galaxies:ISM-galaxies:starburst- ISM:abundance-radio 
hnes:ISM 



1. INTRODUCTION 

The physical properties of molecular clouds are impor- 
tant parameters which are believed to influence largely 
the rate at which stars can form by gravitational col- 
lapse of cloud fragments. The nearby starburst galaxy 
M82 is known to host a large amount of molecular gas 
(M(H2) PS 1 X IO^Mq, e.g. Wild etal. (1992)) and is 
currently undergoing a phase of intense star formation. 
Therefore it has been a prime target to determine the phys- 
ical conditions of the molecular phase of its ISM in order 
to increase our understanding of the relevant physical pro- 
cesses connected with violent star formation activity. 
Most of these studies use CO emission lines to derive the 
physical properties of the molecular gas (e.g. Harris et al. 
(1991); Wild etal. (1992), and more recently Mao etal. 
(2000); Petitpas & Wilson (2000); Weifi etal. (2001)). Yet 
the interpretation of CO emission lines requires the ap- 
plication of radiative transfer models which use simpli- 
fying assumptions about the structure and kinematics of 
molecular clouds. Therefore it is desirable to observe other 
molecules which directly trace the physical conditions and 
thus allow a test of conclusions gained from CO observa- 
tions. Ammonia (NH3) is such a molecule: the relative 
population of its meta-stable inversion lines is sensitive to 
the kinetic gas temperature (e.g. Walmsley & Ungerechts 
(1983)). Therefore observations of ammonia inversion lines 
can be used as a molecular cloud thermometer. Here we 
report the first detection of ammonia towards the most 



prominent molecular source in M82, the south-western 
molecular lobe. We compare our interpretation of the am- 
monia lines to conclusions deduced from CO observations. 

2. OBSERVATIONS 

The (J,K) = (1,1), (2,2), (3,3) and (4,4) inversion lines of 
ammonia were observed towards the south-western (SW) 
molecular lobe in M82 (aaooo = 09''55™49.06'', Jzooo = 
69°40'41".2) using the Effelsbcrg 100m telescope of the 
MPIfR equipped with a dual channel K-band HEMT re- 
ceiver. The system temperature was 180-200 K on a main 
beam brightness temperature (Tmb) scale. The observa- 
tions were carried out in May 2000 and Feb. 2001. The 
beam size of the 100-m telescope at 23 GHz is w 40". The 
data were recorded using an autocorrelator with 8 x 512 
channels and a bandwidth of 40 MHz for each backend, 
leading to a channel spacing of 1 kms~^. The eight back- 
ends were configured such that all four ammonia lines were 
observed simultaneously. Each pair of backends centered 
on a given ammonia line was used to sample both linear 
polarizations independently. To avoid a contamination of 
the data by low-level baseline instabilities we shifted the 
center velocity of the backends (May 2000) and permuted 
the backends (Feb. 2001) every 90 minutes. 
The measurements were carried out in dual-beam switch- 
ing mode, with a switching frequency of IHz and a beam 
throw of 2' in azimuth. Flux calibration was obtained by 
observing both the continuum and the ammonia lines of 
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W3(0H) before each observing run (for fluxes see Mauers- 
berger, Wilson, & Henkel (1988) and Ott etal. (1994)). 
Pointing was checked every 1.5 hours on the nearby con- 
tinuum source 0836+71 and was found to be stable to 
within 5 — 10". A linear baseline was removed from each 
spectrum and intensities were converted to a Tmb scale. 
The summed spectra were smoothed to a velocity resolu- 
tion of 16 kms~^. We estimate the flux calibration of the 
final reduced spectra to be accurate within ±20% (10% 
error of the flux calibration and 10% uncertainty due to 
low-level baseline instabilities). 

3. RESULTS 

The observed spectra of the ammonia lines arc shown 
in Fig. 1. The (J,K) = (l,l)-(3,3) lines are detected with 
a S/N ratio better than 4. The (4,4) line is not de- 
tected. The (1,1) and (2,2) lines peak at visr = 100 kms~ 
(Cioo) and show a weaker component at visr — 160 kms"""^ 
(Ci6o)- Both velocity components are also detected in low 
and mid-J CO emission lines (see e.g. Harris etal. (1991); 
Wild etal. (1992); Mao etal. (2000); Petitpas & Wilson 
(2000)). At 40" resolution the intensity of the low J CO 
transitions from Cieo is stronger than that from Cioo- 
To emphasize this point we show the line profile of the 
^^CO(J =1^0) transition at 40" resolution in the top 
panel of Fig. 1. 

Beam averaged column densities for individual inversion 
states were calculated using 
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(e.g. Henkel etal. (2000)). The column density A^, the 
frequency v and the integrated line intensity are in units 
of cm~^, GHz and K kms~^, respectively. Line parame- 
ters and column densities are summarized in Tab. 1. Note 
that this approximation assumes optically thin emission 
and the contribution of the 2.7 K background to be neg- 
ligible (Tex ^ 2.7 K). Following the analysis described 
by Henkel etal. ((2000)) the rotation temperature {Trot) 
between levels J and J' can be determined from the slope 
of a linear fit in the rotation diagram (normalized column 
density vs. energy above the ground state expressed in 
E/k) by 

a a 
where a is the slope of the linear fit. 

The rotation diagram for the observed ammonia lines 
is shown in Fig. 2. The rotation temperature between 
the (1,1) and (2,2) inversion levels of para-ammonia is 
Tyot = 29K (thick solid line). The uncertainty, derived 
from the extrema of the slope including the (4,4) level as 
an upper limit, is ±5K (dotted lines). The (J,K) = (3,3) 
line is not included in the fit because it belongs to ortho- 
ammonia, thus to a different ammonia species. Neverthe- 
less, it nicely fits to the rotation temperature derived from 
the two lowest inversion levels of para-ammonia. 
In addition to the analysis of the integrated line intensities 
we have decomposed the (J,K) = (1,1) and (2,2) spectra 
into two Gaussian components with fixed center veloci- 
ties of Visr = 100 kms~^ (Cioo) and 160 kms~^(Ci6o)- 
The Gaussian decomposition is shown together with the 



observed line profile in the two top panels of Fig. 1. A sep- 
arate analysis of the rotation temperatures between the 
(1,1) and (2,2) inversion levels for both components yields 
Trot, 100 = 24tg^,K and TrotMO = 'HtHK for Cioo and 
Cioo, respectively. For the errors we have assumed 30% 
uncertainty from the Gaussian decomposition. Note that 
the observed upper limit for the (J,K) = (4,4) line inten- 
sity is still consistent with the shallow slope derived for the 
emission arising from the 160 kms^^ component (Cioo)- 

4. DISCUSSION 
4.1. The NH3 emitting volume 

High-spatial resolution CO observations show that the 

emission at visr ~ 100 kms~^ is mainly associated with 
the SW molecular lobe. The emission at visr ~ 160 kms~^ 
arises from regions closer to the nucleus , i.e. from the 
central molecular peak and the inner CO outflow region. 
Thus C160 covers the regions of the western mid infrared 
peak close to the central molecular peak (Tclcsco & Gezari 
(1992)), and the region where an expanding molecular su- 
perbubble was identified (Weifi etal. (1999)). Therefore 
C160 represents a more active region than the outer molec- 
ular lobe itself. An overlay of the two regions at high spa- 
tial resolution (w 2") as observed in the ^^C0(J = 2— >1) 
transition (Weifi etal. (2001)) is shown in Fig. 3. 

4.2. Comparison with temperature determinations from 
other line observations 

Radiative transfer calculations (e.g. Walmsley & Un- 
gerechts (1983)) of NH3 show that rotation temperatures 
determined from meta-stable levels only reflect the kinetic 
gas temperature for low (T < 15 K) temperatures. For 
larger temperatures the rotation temperature largely un- 
derestimates the kinetic gas temperature due to depopula- 
tion mechanisms. Correcting our results for these effects, 
our mean rotation temperature of 29 K corresponds to a 
mean kinetic temperature of Tkm ~ 60 K. Using the Gaus- 
sian decomposition we find Tfej„ « 45K for the tempera- 
ture in the SW molecular lobe (Cioo), and Tkm ~ 80K for 
the regions closer to the nucleus (Cieo). 
These values are in good agreement with recent high- 
resolution kinetic temperature estimates using radiative 
transfer calculations based on CO observations by Weifi 
et al. ((2001)). They derive Tki„ « 55 K for the SW molec- 
ular lobe, and an average of Tkin ~ llOK over the region 
which corresponds to Cieo m our study. The simultane- 
ous observations of two [C I] fine structure lines towards 
the SW lobe allowed Stutzki etal. ((1997)) to derive a 
lower limit for the kinetic temperature of the [C I] emit- 
ting gas in Cioo of 50 K which is consistent with the tem- 
perature derived above from the NH3 lines. Similar val- 
ues {Tkin = 50 K) have been found by Seaquist & Frayer 
((2000)) using HCO+ and HON emission lines. 
Note, however, that the uncertainty of the kinetic temper- 
ature estimates from the Gaussian decomposition of the 
NH3 spectra is quite large. Nevertheless, the results are 
in line with the general picture that the kinetic tempera- 
ture rises from the outer parts of the SW molecular lobe 
towards the active regions closer to the nucleus of M82. 
In comparison with rotation temperatures determined 
from NH3 in IC342 {Trot « 50 K; Martin & Ho (1986)) 
and Maffei 2 {Trot ~ 85 K; Henkel etal. (2000)) the values 
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determined here for M 82 are very low. Thus the tempera- 
ture of the dense moleeular gas in M 82 is obviously much 
lower than in other galaxies with large nuclear concentra- 
tions of molecular gas. This explains the apparent un- 
derabundance of molecules like SiO, CH3OH, HNCO, and 
CH3CN in M82 which only form in a dense and warmer 
environment (see e.g. Mauersberger & Henkel (1993)). 

4.3. NH3 abundance 

Using eq. A15 of Ungerechts, Walmsley, & Winnewisser 
((1986)) we have estimated a total beam-averaged NH3 
column density of iV(NH3) w 1 x lO^^cm"^. The beam- 
averaged H2 column density was estimated to be A^(H2) « 
1.9 X 10^^ cm^^ using the CO spectrum shown in Fig. 1 and 
a conversion factor of Xqo = 5 x 10^^ cm~^ (Kkms~^)~^, 
which corresponds to the average value of the conversion 
factors derived by Weifi et al. ((2001)) in the region covered 
by 100 180 kms""'^ emission. This yields a relative abun- 
dance of ammonia of X(NH3) « 5 x 10^^", which is an 
extremely low value. In nearby dark clouds the fractional 
NH3 abundance is of order X(NH3) ss lO"'^ (e.g. Benson 
& Myers (1983)) and in hot cores NH3 is even more abun- 
dant (X(NH3) w lO"-"^- -"^; e.g. Mauersberger, Henkel & 
Wilson (1987)). In a recent work on Maffei 2 Henkel etal. 
((2000)) found XCNRs) « IQ-^. 

A closer look at the (J,K) = (1,1) and the ^^COiJ = 
1^0) spectra displayed in Fig. 1 (top) reveals an- 
other interesting aspect regarding the spatial variation 
of the NH3 abundance: the line temperature ratio 
T(NH3(l,l))/T(i2CO(J=1^0)) decreases with increas- 
ing velocity for visr > 100 km s"""^. The solid-body ro- 
tation in the inner part of M82 allows one to asso- 
ciate velocities lower than the systemic velocity of visr = 
230kms~^ with a specific distance from the nucleus (e.g. 
Neininger etal. (1998)). Therefore, T(NH3(1,1)) per ve- 
locity interval traces the ammonia column density in a 



region much smaller than the spatial resolution, while 
T(i2CO(J=1^0)) traces the H2 column density. Note 
that conversion from I(CO) to iV(H2) changes across the 
major axis of M82 and therefore for each velocity inter- 
val (Weifi etal. (2001)). Fig. 4 shows a histogram with 
Ar(NH3) versus radial velocity and galactocentric radius, 
accounting for variations of Xqo as determined by Weifi 
etal. ((2001)). The NH3 abundance is found to decrease 
towards the center of M82. We believe that this find- 
ing reflects a real change of the ammonia abundance, be- 
cause a constant NH3 abundance would imply that Xqo 
changes by a factor of more than 40 on a linear scale of 
300 pc which is not consistent with radiative transfer mod- 
els. Due to its low energy threshold for photodissociation 
(« 4.1 eV, Suto & Lee (1983)), NH3 should be destroyed 
rapidly in PDRs (Glisten & Fiebig (1988)). This process 
should be even more efficient when the bulk of the gas is 
distributed in a diffuse phase with low H2 column densities 
as it seems to be true for the central regions in M 82 (Mao 
etal. (2000), Weifi etal. (2001)). In such an environment 
shielding against the strong UV radiation is inefi:ective 
which leads to low NH3 abundances towards the central 
star forming regions. Note, that NH3 abundances are not 
always small in a warm environment. Efficient release of 
NH3 into the gas phase by evaporation of dust grain man- 
tles is known to occur in galactic 'hot cores'. This process, 
however, becomes efficient at slightly higher temperatures 
than those obtained by us for NH3 in Sect. 4.2. 
We therefore interpret the abundance gradient as a result 
of ammonia being almost completely dissociated in the 
harsh environment close to M82's nucleus, whereas it can 
still partially survive, with low abundance, in the cooler, 
denser, and thus better shielded, clouds of the SW lobe. 

This research project was supported by the Deutsche 
Forschungsgemeinschaft grant SFB 494. 
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Fig. 1.— NH3 spectra towards the SW part of M82 (02000 = 09''55'"49.06^ (520oo 
for each spectrum. The high-resolution spectrum in the top panel shows the ^^CO(J = 
resolution (40") scaled by 1/370. For the fits to the (1,1) and (2,2) lines see Sect. 3. 



3°40'41".2). The velocity resolution is 16 kms"^ 
> 0) emission line at the same position and spatial 
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Fig. 2. — Rotation diagram of meta-stable ammonia transitions towards the SW molecular lobe in M82. The filled squares show the 
normalized column densities determined from the integrated line intensities. The thick solid line (denoted by 'total') corresponds to a linear 
fit to the (J,K) = (1,1) and (2,2) lines for these values. The dotted lines correspond to the linear fits with the lowest and highest slope, which 
are still consistent with the data including the upper limit for the (J,K) = (4,4) line. The open squares and the stars show the normalized 
column densities determined from a Gaussian decomposition of the (J,K) = (1,1) and (2,2) spectra into the velocity components Cioo and 
Ci60- Note that for display purpose the values for Cieo and Cioo have been shifted by -0.2 and -0.8 on the y-axis scale, respectively. The 
thin solid lines correspond to linear fits to these data points. 
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Fig. 4. — NH3 abundance relative to H2 per velocity interval. The upper axis denotes the distance of the emission region from the nucleus 
of M82 assuming pure solid body rotation and a velocity gradient of 0.7 kms~^ pc~^. 
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Table 1 

Parameters of the ammonia lines towards the SW molecular lobe of M82. Column 2 to 6 correspond to the 

LINE PARAMETERS OF THE ENTIRE LINE PROFILES DERIVED FROM A VISUAL INSPECTION. COLUMNS 7 AND 8 GIVE THE COLUMN 
DENSITIES PGR THE (1,1) AND (2,2) TRANSITIONS DERIVED FROM A GAUSSIAN DECOMPOSITION OF THE SPECTRA. 



Transition 

(J,K) 


rppeak 
-'-mb 

(mK) 


/ Trabdv 

(K kms-i) 


T 7-pea/c 
^LSR 

(kms^i) 


m/2 

( km s 1 ) 


N{J, K) 
(10^2 cm-2) 


N{J,K)ioo 
(10^2 cm-2) 


iV(J,i^)l60 

(10^2 cm-2) 


(1,1) 


6.9 ±0.8 


0.61 ±0.16 


97± 16 


94± 11 


4.0 ± 1.1 


2.2 ±0.8 


1.8 ±0.6 


(2,2) 


4.3 ±0.9 


0.36 ±0.11 


105 ± 16 


80± 10 


1.8 ±0.5 


0.9 ±0.3 


1.0 ±0.3 


(3,3) 


3.5 ±0.7 


0.14 ±0.06 


95 ±10 


47 ±20 


0.6 ±0.26 






(4,4) 


< 2.1(3(7) 


< 0.1 






< 0.4 







